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Abstract 

Magnetohydrodynamics of strongly magnetized relativistic fluids is derived 
in the ideal and dissipative cases, taking into account the breaking of spa- 
tial symmetries by a quantizing magnetic field. A complete set of trans- 
port coefficients, consistent with the Curie and Onsager principles, is de- 
rived for thermal conduction, as well as shear and bulk viscosities. It is 
shown that in the most general case the dissipative function contains five 
shear viscosities, two bulk viscosities, and three thermal conductivity coef- 
ficients. We use Zubarev's non-equilibrium statistical operator method to 
relate these transport coefficients to correlation functions of the equilibrium 
theory. The desired relations emerge at linear order in the expansion of the 
non-equilibrium statistical operator with respect to the gradients of relevant 
statistical parameters (temperature, chemical potential, and velocity.) The 
transport coefficients are cast in a form that can be conveniently computed 
using equilibrium (imaginary-time) infrared Green's functions defined with 
respect to the equilibrium statistical operator. 
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1. Introduction 

Computation of transport coefficients of relativistic systems is of impor- 
tance in a variety of fields, including relativistic astrophysics, cosmology, and 
heavy-ion physics. The kinetic theory based on the Boltzmann equation pro- 
vides an adequate tool for describing transport in these systems in real time 
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in the dilute regime. It may be systematically extended to denser systems 
in cases where one can establish a hierarchy of the timescales governing the 
system. An alternative to the Boltzmann approach is based on the Kubo 
method, in which non-equilibrium processes are regarded as the response 
of the system to an external perturbation, which can be computed from 
equilibrium correlation functions. This method has the advantage that the 
transport coefficients can be evaluated by using equilibrium Green's functions 
formulated in imaginary time. 

In this paper we shall adopt the correlation-function method developed by 
Zubarev [1, 2, 3] to derive expressions for transport coefficients of strongly 
magnetized relativistic fluids. The Zubarev method is based on the con- 
cept of a non-equilibrium statistical operator, which is a generalization of 
the equilibrium Gibbs statistical operator to non-equilibrium states. This 
approach will enable us to obtain kinetic transport coefficients in terms of 
correlation functions by expanding the operator in terms of small gradients 
of the thermodynamical parameters. It is generally assumed that the state 
of the system can be defined in terms of (space-time dependent) fields of 
temperature, chemical potential, and momentum — a typical parameter set 
that is sufficient for a complete description of the system. (The set of relevant 
parameters depends in general on the specific problem under consideration.) 
The transport coefficients are then obtained from linear perturbations of the 
non-equilibrium statistical operator around its equilibrium value. In this pa- 
per we will be concerned with the situation in which the system is in an 
external Abelian gauge field, i.e., an electromagnetic field. 

The physical motivation for studying relativistic fluids in strong magnetic 
flelds is provided by compact stellar objects (neutron stars). It is now well 
established that neutron stars feature magnetic fields of order 10^^-10^^ G 
on average, with a subclass of neutron stars (magnetars) having magnetic 
fields of order 10^^-10^^ G. Although the origin of such large magnetic fields 
is not known, one possibility is the conservation of magnetic flux [4, 5] under 
collapse from ordinary stellar dimensions. Alternatively, the field could be 
generated by a dynamo mechanism [6]. The field strengths B < 10^^ G 
inferred for magnetars [7, 8] refer to the surface of a compact star. The 
internal fields could be by several orders of magnitude larger. An upper limit 
on the magnitude of the magnetic field B < 10^^ G is set by the virial theorem 
for a star in gravitational equilibrium [9, 10]. Observationally, soft gamma- 
ray repeaters (SGRs) and anomalous x-ray pulsars (AXPs) are now identified 
with magnetars; these objects have been associated with specific supernova 
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remnants and form a subclass of young strongly magnetized neutron stars 
which are distinct from ordinary radio pulsars. 

The equilibrium properties of dense matter in strong magnetic fields have 
been studied intensively since the discovery of magnetars. These properties 
include the equation of state of dense nuclear and quark matter [11, 12, 13], 
neutrino propagation and pulsar kicks, and neutrino emission and cooling [14, 
15, 16, 17, 18, 19, 20, 21]. However, hydrodynamics and non-equilibrium 
processes in strong fields have not been studied until recently. 

In an earlier paper [22], we have developed an anisotropic hydrodynamic 
theory to describe strongly magnetized strange matter in a compact star. It 
was shown that the dissipative processes in matter are completely described 
in terms of seven viscosity coefficients (which include five shear viscosities 
and two bulk viscosities), three thermal conductivities, and three coefficients 
of electrical conductivity. Our model is the relativistic counterpart of the 
non-relativistic formulation of plasma theory in strong magnetic fields due 
to Braginskii [23]. 

The purpose of this work is to extend the description of the hydrodynam- 
ics of strongly magnetized fiuids by relating the hydrodynamic equations to 
certain correlation functions. The latter can be computed using the meth- 
ods of equilibrium statistical mechanics. Among several methods to derive 
Kubo formulas, we choose the one based on the idea of non-equilibrium sta- 
tistical operators (sec Refs. [1, 2, 3] for a detailed discussion of this method 
for non-magnetized fluids.) Thus, compared to other relativistic treatments 
of magnetohydrodynamics (see, e. g., Refs. [24]), our approach relates the 
transport coefficients entering magnetohydrodynamics to retarded correla- 
tion functions of the microscopic theory. 

This paper is organized as follows. Section 2 is devoted to magnetohy- 
drodynamics of relativistic fluids in the ideal (Subsec. 2.1) and dissipative 
(Subsec. 2.2) cases. We derive a thermodynamically consistent set of equa- 
tions describing fluid dynamics in strong magnetic fields and identify the full 
set of transport coefficients needed for a complete description of dissipation. 
In Sec. 3 we utilize the non-equilibrium statistical-operator method in the 
formulation by Zubarev to derive Kubo formulas relating the transport coef- 
ficients appearing in our magnetohydro dynamic theory to correlation func- 
tions defined in terms of the underlying elementary fields. A discussion of 
our results and a comparison to previous work is given in Sec. 4. Natural 
units h = kB — c — 1 are adopted, and the SI system of units is em- 
ployed in equations involving electromagnetism. The metric tensor is given 
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by 5'^^ = diag(l, —1, —1, —1). A number of calculational details are relegated 
to the appendices. 

2. Magnetohydrodynamics of relativistic fluids 

2.1. Non-dissipative fluids 

Hydrodynamics arises as a long- wave-length, low-frcqucncy effective the- 
ory of matter from an expansion of the energy-momentum tensor T^^", the 
current (s) n'^ of conserved charge(s), and the entropy density flux s'^, etc., 
with respect to (small) gradients of the four-velocity u^^ and the thermody- 
namic parameters of the system, such as the temperature T, the chemical 
potential /j, etc. The zeroth-order terms in this expansion correspond to an 
ideal fluid, which we will discuss here; the extension to the dissipative case 
is given in the following subsection. 

For the sake of simplicity, we will consider only one conserved charge, 
e. g., the electric charge. Many physical systems of interest may have other 
conserved charges, such as baryon number, isospin, etc. The extension of 
our discussion to multiple conserved charges is straightforward. The hydro- 
dynamic equations can be expressed as conservation laws for the energy- 
momentum tensor T^^'^ and the conserved currents, here the electric current 



d,T^- = 0, (1) 
= 0. (2) 

In the presence of an electromagnetic field the zeroth-order contribution to 
the energy-momentum tensor can be written as [22, 25, 26] 

-'o ~ ^Fo 1" -'em? K"^) 
T^^ = £M^M^-PA^^-i(M'^^F;,^ + M^^F^^), (4) 

TZa = -F^^F\ + ^—F'"^F,„, (5) 

where F^^ is the field- strength tensor, e and P are the local energy density 
and the thermodynamic pressure, A^*^ = g^^'^ — u^u^ is the projector on the 
directions orthogonal to m^, and M^'' is the polarization tensor. In the non- 
dissipative limit, the entropy is conserved in addition to the charge, i.e., it 
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obeys a continuity equation analogous to Eq. (2). The charge and entropy 
currents in the non-dissipative theory are 

< = nu^, (6) 
s^o = su^, (7) 

where n and s are the electric charge and entropy densities measured in the 
rest frame of the fluid. For the following discussion it will be convenient to 
decompose the tensor F^'^ into components parallel and perpendicular to u^, 

= E^'u" -E''u^ + ^e''''^l'{uc.B^-upB^), (8) 

where we have defined the four-vectors E^^ = F^^u^ and = e^^'^^ F^a^p/2, 
with e^"°'^ denoting the totally anti-symmetric Levi-Civita tensor. In the 
rest frame of the fluid, li'* = (1, 0), we find that E^ ^ B^ ^ 0, E' = and 
B^ — —e^^^Fjk/2, where the indices i, j, /c run over 1, 2, 3. The latter three 
vectors are precisely the electric and magnetic fields in this frame. Therefore, 
E^ and B^ can be interpreted as the electric and magnetic fields measured 
in the frame in which the fluid moves with a velocity u^. 

The antisymmetric polarization tensor M^'' describes the response of mat- 
ter to the applied field strength F^^ . In the case of systems described by a 
grand canonical ensemble, it is given by M^^ = —dQ/dFn^ in terms of the 
thermodynamic potential ^2. For later use, we also define the in-medium field 
strength tensor H^^'^ = F'^^ — M^^. In analogy to what is done for F''^, we 
can decompose M^^ and H'^'^ as 

Mi^y = P^'u" - P^^u" + ]^e^''"'^ {Mpucc- M^Ufi), (9) 

i/M- ^ Di^u- _ D'^u'^ + L'^'"''^ {H^Ua- H^u^) , (10) 

with P^ = -M^"'u^, M*" = e'"'"^^ M^aU0/2, H^' = e'"'"^^ and Di" = 

In the rest frame of the fluid, the non-trivial components of these ten- 
sors are (F^^, F^o, F^O) = E, (F^^, F^^, F^i) = B, (M^o, M^o, M^O) = -P, 
(M32,Mi3,M2i) = M, (//lo, ^30) ^ and {H^\H^\H''^) = H. 
Here P and M are the electric polarization vector and magnetization vec- 
tor, respectively. In linear approximation they are related to the flelds E 
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and B by the usual expressions P — XeE and M = x^B, where Xe ^-nd 
Xm are the electric and magnetic susceptibilities. Note that the four-vectors 
E^, ^ ■ ■ ■ are all space-like, i.e., E^u^ = 0, B^u^ = 0, ■ ■ ■, and are normal- 
ized as E^E^ = -^2, B^'B|, = -S^, • • •, with E = |E| and B = \B\. 

For many systems encountered in astrophysics, one example being the 
interior of a neutron star, the electric field is much weaker than the magnetic 
field. Therefore, in the following discussion we will omit the contributions 
originating from the electric field, but will occasionally comment on how 
our expression will be modified by it. Next let us introduce the four- vector 
b'^ = B^/B, which is normalized as h^hfj, = — 1 (in agreement with the nor- 
malization of the vector B^^ above), along with the antisymmetric rank-2 
tensor h^^ = e^^^^h^up. Then, upon dropping the electric field in the field 
strength tensors, we write these in the compact form 

F^'^ = -m^^ (11) 
Ml''' = -Mfc'^^ (12) 
H'''' = -Hh''\ (13) 

where M = |M| and H = |H|. Again neglecting electric fields, the mat- 
ter and field contributions to the energy-momentum tensor (3) can now be 
written in terms of b^^ and b'^'^ as (see, e.g., Refs. [22, 24]) 

T^^ = su^'u'' - P^E"" + Piib^b" , (14) 

T^M = iB^uV-Ei^'^-b^b''), (15) 

where S'''^ = A'^'^ + b'^b'^ is the projection tensor on the direction orthogonal 
to both u^^ and 6^. Here we have defined the transverse and longitudinal 
pressures P± — P — MB and P\\ — P relative to the vector b^^. In the 
absence of a magnetic field, the fiuid is isotropic and P± — P\\ — P, where 
P is the thermodynamic pressure defined in Eq. (5). In the local rest frame 
of the fluid, we have 6^ = (0,0,0,1), where without loss of generality the 
direction of the magnetic field is chosen as the 2;-axis. We see then that the 
electromagnetic tensor takes the usual form, while Tp^ = diag(e, P^, Pi, Py), 
as required. 

Next we would like to check the consistency of the terms that appear in 
TpQ with the formulas of standard thermodynamics involving electromagnetic 
fields. By using the thermodynamic relation 

e = Ts + nn-P (16) 
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and the conservation equations for tIq and Sq i'^ ideal hydrodynamics, one 
can show that the hydrodynamic equation Uj^d^T^^ — together with the 
MaxweU equation (22) imphes 

De = TDs + ^lDn - MDB, (17) 

which is consistent with the standard thermodynamic relation 

de = Tds + ndn- MdB. (18) 

Prom Eq. (16) and Eq. (18), one can obtain the Gibbs-Duhem relation, 

dP = sdT + ndn + MdB. (19) 

We should note that the potential energy —MB has already been included 
in our definition ole. Otherwise, new terms —MB, —D{MB), and —d{MB) 
should be added to the left-hand sides of Eq. (16), Eq. (17), and Eq. (18), 
respectively, while the Gibbs-Duhem relation keeps the form of Eq. (19). 

The complete set of non-dissipative hydrodynamic equation includes the 
Maxwell equations, which we state for completeness: 

d,{BV - B''u>') = 0, (20) 
df.H^'' = n\ (21) 

Gontracting Eq. (20) with gives the so-called induction equation 

e + DhiB- u^lfd^K = 0, (22) 

where 9 = d^v/^ is the velocity divergence and D = u^d^j^ is the substantial 
(co-moving) time-derivative. (Note that b^Db^ = b^^Dbf^ = 0.) In closing this 
section, we confirm that the form of the energy -momentum tensor TpQ is con- 
sistent with standard thermodynamic formulas for matter in electromagnetic 
fields [22]. 

2.2. Dissipative fluids 

The next-to-leading (first-order) derivative expansion of conserved quan- 
tities leads to the Navier-Stokes-Fourier-Ohm theory which, as is well known, 
includes dissipative effects. The continuity equation for the entropy now 
changes to 

Td^s" > (23) 
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(second law of thermodynamics). In the dissipative theory the quantities 
j'nv ^ n'', and can be generally expressed as 

= su^ + f,, (24) 

where /i'^ is the energy flux, r''^ is the viscous stress tensor, and j'^ and 
arc the charge and entropy diffusion fluxes. They are all orthogonal to u^, 
reflecting the fact that the dissipation in the fluid should be spatial. We shall 
assume that can be expressed as a linear combination of and j^. This 
allows us to incorporate the fact that the entropy flux is determined by the 
energy and charge diffusion fluxes. Thus, following Refs. [27, 28], we write 

= ih^ - o^f, (25) 

the coefficients 7 and a being functions of thermodynamic variables. 

In order to discuss the dissipative parts, let us first define the four- velocity 
u'^, since it is not unique when one allows for energy exchange by thermal 
conduction. We employ the Landau-Lifshitz frame, in which u'^ is chosen 
to be parallel to the energy density flow, i.e., u^T'^" = u^e. It follows from 
the first equation (24) that = 0. We project Eq. (1) onto u'^ and exploit 
the fact that electric field Ex is zero. Utilizing Eq. (25), straightforward 
manipulations lead to 

Ds+{£ + P)e- T^^d^u^ + MDB = j^Ex = 0. (26) 

Combining this equation with the thermodynamic relation e = Ts + iin — P 
and the continuity equation (2), we arrive at 

Td^s'' = T''''w^,~rTV.a + {i2-Ta)d^r, (27) 

where Vo- = A^^d'^ and uf^ = | (V'^u'^ + V^u"). The first and second terms 
on the right-hand side of Eq. (27) (the dissipative function) correspond to 
viscous and thermal dissipation, respectively, whereas the last term vanishes 
due to our choice of a. For a thermodynamically and hydrodynamically 
stable system the dissipative function should be non-negative. This implies 
that for small perturbations it must be a quadratic form, so we can write 

Ta = /X, (28) 
= 7/^^"^m;«;3, (29) 
f = K^'^TW^a, (30) 
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where r)'^'^"'^ is the rank-four tensor of viscosity coefficients and n^'^ is the ther- 
mal conductivity tensor with respect to the diffusion flux of electric charge 
density. By definition, r^'^'^"^ is symmetric in the pairs of indices a, (5 and 
/i, I/. It necessarily satisfies the condition fjt^""!^ (^B'^) = ri'^'^^^^—B"), which is 
Onsager's symmetry principle for transport coefficients. Similarly, the ten- 
sors and a^^^ should satisfy the conditions k^^{B^) — k^'^{—B^) and 
a^'^{B^) = a'^^{—B^). Furthermore, all the tensors of transport coefficients 
^fiua/B^ i^fiu ^ ^fiu ]^^-^g^ orthogonal to -u^ by definition. 

The constructive equations connecting the irreversible fiuxes r^,y , j^u and 
the thermodynamic forces W/j^i, , VjuCc in dissipative magnetohydro dynamics 
may have a more general tensor structure 

f = ^'''wp. + K^^TV^a, (31) 

where X^^p and 7^^"^ are rank-three tensors representing the coupling between 
charge diffusion and energy- momentum transport. However, the Curie prin- 
ciple, which states that linear transport relations can be realized between 
irreducible tensors of the same rank and parity, does not allow for such higher- 
rank tensors X'^^'^ and 7'*''^. There are, however, some exceptions as, e. g., 
systems with parity violation [29]. In strong magnetic fields the ranks of 
the tensors are changed compared to the non-magnetic case because of the 
breaking of spatial symmetries by the magnetic field. However, the parity 
of the tensors with respect to reflections should be preserved. We therefore 
conclude that X^'^p = = ■ 

The tensors rj^^'^l^ and k^^ can be decomposed into sums of irreducible 
tensor combinations constructed from the vectors m^, 6^, g^'^ ^ and 6^*". In the 
case of the tensor rj^''^'^^, all the independent irreducible tensor combinations 
that share its symmetries and are orthogonal to u^' are given by [22] 

(i) A^^A""^, 

(ii) A^^A'^^ + A^'^A''", 
(in) A^'^b'^b'^ + A"^6'^6^ 

(iv) bf'b^b^b'^, 

(v) A^^^bi^ + A'^^b^b'' + A'^^fo^fe" + A'^W^, 

(vi) A''"^^^ + A^^6^° + A'^^fe^" + A'^'^b'''^, 

(vn) 6'^" 6^6^ + b'^^f'b'' + bP^b^b"" + (32) 
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In Ref. [22] the combination h'^°'}f^ + h'^l^lf°' was added to the above hst. 
However, as shown in Appendix Appendix A, this tensor can be written as 
a hnear combination of tensors hsted in Eq. (32); specificaUy h^^h^P -\-l)^^l^}f°' — 



The independent irreducible tensor combinations that have the symmetry 
of k'^^ and are orthogonal to read 



In order to separate the contributions of shear and bulk viscosities we will 
need to construct tensor structures that differ from those used in Ref. [22]. 
An appropriate set is given by 



(i') Ai'^A'^P = (i), 

(ii') Af^A"^ + Af'^A'"^ = (ii) , 

(iii') S'^'^S"'' = 2 (iii) + 2 (i) + 2 (iv) , 

(iv') IfU^h'^h^ = (iv), 

(v') 'E^^U'h^ + 'S'^h^'h'^ + 'E^'^U'h'^ + El^^lfy^ = (v) + 4(iv) , 

(vi') E'^lf^ + Ey^h^^ + B'^lf'^ + Ey^h^^ = (vi) + (vii) , 

(vii') h'^U'hf^ + U'^hi'h'' + hi'^U'h'' + U^'^h^'l/ = (vii) . (34) 



Prom these tensors we construct the following combinations, which are either 
traceless. 



2(i)-(ii) + 2(iii)-(v). 



(viii) A' 



(ix) 6^6^ (x) 



(33) 



(i") A'^^A'^^ + A^'^A^ 



l^lJ-v 




(iii") - E^'^Vh^ - E^'^^b" - E'^^U'h'^ - B 



(v") b^'^b'^b'^ + b^'^lfb'^ + If^b^'b'^ + b'^'^lfb^, 



(35) 



or have non- vanishing trace. 



(vi") ^S'^'^S'^^ 
(vii") Zb^b'^b'^b^ . 



(36) 
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Now we can write down the most general expression for r)^^'^^ as 

Similarly, we also have the general decompositions for k'^" , namely 

k"" = K^E"'' - K\ilfb'' - K^b"'' . (38) 

Thus we conclude that a fluid in a magnetic flcld has in general seven in- 
dependent viscosity coefficients and three independent thermal conduction 
coefficients, which are defined through Eq. (37) and Eq. (38), respectively. 
Substituting Eqs. (37) and (38) into Eqs. (29) and (30), respectively, we 
obtain the expressions for the viscous stress tensor and thermal flux in terms 
of viscosity coefficients and thermal conductivities, 

- 2r]2 {If^y^h^ + U'Ei^h^) - 2773 (S^"&'^^ + Ey^hi"^) w^p 

+ 27?4 [ir'h^^ + h'^^h^^hf") + ^CxS'^"'/' + 3C||&^fe> , (39) 



f = K^TEi'^V'^a- K^b^b^TV^a- K^b'^'TVua , (40) 

where = S'^'^w^j^ and (p = b^b^w^i^. The coefficients of the tensor r'^'^ 
multiplying traceless tensors are identifled as shear viscosities, while the co- 
efficients multiplying the tensors with non-vanishing trace are identifled as 
bulk viscosities. The k's are the thermal conductivities. 

The divergence of the entropy-density flux (27) can now be written ex- 
plicitly as 

+ 2772 {b^bX" - b'^bXn ip^bPwp, - b^lfwp^) 
+ \Ci.4>^ + K\\^^-^^iT^E>^^V,aV,a + K\\T^{lfV,af. (41) 

One should note that the terms corresponding to the transport coefficients 
r]3,r)4, and in Eqs. (39)-(40) do not contribute to the divergence of the 
entropy-density flux. For stable systems, all the other transport coefficients 
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must be positive to satisfy the second law of thermodynamics (positivity of 
the dissipative function.) There are, however, examples in which this condi- 
tion may be violated. Indeed, as shown in Ref . [22] , strange quark matter (if 
it occurs in some form in compact stars) may have negative transverse bulk 
viscosity C_l in strong magnetic fields, in situations where non-leptonic weak 
processes are dominant. The change in the sign of this transport coefficient 
can be linked to the onset of a hydrodynamic instability [22]. 

Having determined the general form of the relativistic hydrodynamic 
equations in strong magnetic fields, we turn now to the problem of determin- 
ing the transport coefficients entering these equations. We have previously 
determined the bulk viscosity [22] in the specific case of strange matter, where 
knowledge of the equation of state and pcrtTu4:)ative weak interactions rates is 
sufficient to fix the value of (. In order to obtain the remaining transport co- 
efficients one must resort to the microscopic theory. Two general approaches 
are available for this task: one is based on the solution of a (hnearized) Boltz- 
mann equation; the second is based on the evaluation of certain correlation 
functions. The first approach is suitable for dilute systems or systems with 
well-defined quasiparticles. The second approach is more general, since the 
required correlation functions can in principle be evaluated for a system with 
arbitrarily strong interactions far from equilibrium. We will follow the second 
approach to express the transport coefficients of a fluid in a strong magnetic 
field in terms correlation functions. In order to do so, we apply the powerful 
method of non-equilibrium statistical operators developed by Zubarev [1] . 

3. Non-equilibrium Statistical Operators and Kubo Formulas 

Zubarev's method of non-equilibrium statistical operators (NESO) allows 
one to develop a statistical theory of irreversible processes via a generalization 
of the Gibbs approach to equilibrium systems. Zubarev's method assumes 
that one can construct a statistical ensemble representing the macroscopic 
state of the system in non-equilibrium. In general this is possible if one is 
interested in quantities that are averaged over time intervals that are long 
enough that the initial state of the system is unimportant and the number 
of parameters describing the system is correspondingly reduced. 

We consider the system to be in the hydrodynamic regime (in which the 
thermodynamical parameters such as temperature and chemical potentials 
can be defined locally) and, following Zubarev, choose the NESO in the 
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form [1] 



p{t) = Q ^exp 



Q = Tr exp 



(42) 
(43) 



where the operator Z reads 



Zi^,t) = e[ dtie^(*i-*) [/3"(x,ti)fo,(x,ti)-a(x,ti)7V°(x,ti)] ,(44) 

J —oo 



with e — > +0 after the thermodynamic hmit is taken, and 

r(x,t) = /3{^,t)u%^,t), 
q;(x, t) = /3(x, t)/i(x, t). 



(45) 
(46) 



Physically the parameters ,5, //, and u/* stand for the inverse local equilibrium 
temperature, chemical potential, and flow four-velocity, respectively. The 
conditions needed to make such identifications will be addressed below. The 
operators of the energy-momentum tensor T^^'^ and the electric current TV'* 
satisfy the local conservation laws 



di,T^"' = 0, di,m = 0. 



(47) 



Integrating Eq. (44) by parts and using the conservation laws (47) , we rewrite 
the exponent of Eq. (42) as 



d^xZ(x, t) 



d^x 



/3'^(x,t)To,(x,t)-a(x,t)7V°(x,i) 



f^,(x,ti)5^/3'^(x,ti) 



- N''{K,ti)d^a{K,ti) 



(48) 



The quantities d'^P'^, d^a are thermodynamic "forces" as they involve gradi- 
ents of velocity field, temperature, and chemical potential. It is natural to 
identify the first integral in Eq. (48) as the equihbrium part of the statistical 
operator and the second as the non-equilibrium part of the statistical oper- 
ator (the latter vanishing when the space-time variations of the parameters 
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a{yi,t) and /3'^(x, t) vanish.) Thus we may define the local equilibrium part 
of the statistical operator as 



Pi{t) = Q/cxp 
Qi ^ Tr exp 



rf^x (|/3^(x, t)fo^(x, t) - a(x, t)iV°(x, t)) 
rf^x (^/3'^(x, t)ro,.(x, i) - q;(x, i)7V°(x, i)) 



and write the complete statistical operator as 

Q = Tre-^+^, 

having introduced the following short-hand notation for the integrals 
A{t) - J £x (/3''(x, t)fo,(x, t) - a(x, t)iV°(x, t)) , 



with 



J J —oo 



^(x, t) = T^,(x, t)a'^/3''(x, t) - iV'^(x, t)9^a(x, t). 



(49) 
(50) 

(51) 
(52) 

(53) 

(54) 

(55) 



Now we consider the case of small perturbations from the equilibrium 
state, i.e., the case in which the thermodynamic forces are sufficiently small 
that the time-scales for relaxation of non-equilibrium states are large. In 
this regime, the relations between the thermodynamic forces and irreversible 
currents arc linear, so that we can approximate the non-equilibrium statistical 
operator by keeping only the linear term in the expansion of the exponent in 
Eq. (51). Thus, in the linear approximation the complete statistical operator 
reads 

• 1 



1+ / dr (e-^^Be^^ - {B), 



Pi, 



(56) 



where {B)i = Tr yPiBj is the average over the local equilibrium operator. 

Now we are in a position to evaluate the energy-momentum tensor aver- 
aged over the non-equilibrium distribution. We obtain 



(f^^(x,t)) = Tr p(t)f^'^(x,t) 



(57) 
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where 



Tr 
Tr 



dTe-^^Be'^^pi{t)f'"'{^,t) 



(58) 



is the deviation from the average equihbrium value of the energy-momentum 
tensor (T'^'^(x, t));. Substituting the integral (54) into Eq. (58), we obtain 

J^'dT{f'^''{x,t) [e-^^C'(xi,ii)e^^ - {C{K,,h))i'\)i .(59) 

The notation can be compactified by introducing Kubo correlation functions 
defined as 



X(x,t),f(x',t') 



rfr(X(x,t) e-^-F(x',t')e^^-(F(x',t'))/ )i, (60) 



in terms of which we finally obtain 

T^'^ ^S{f^''{^,t)) = Jd'^ij' dtie^(*i-*)(f^^(x,t),C'(xi,ti)).(61) 

Similarly, the linear response of the charge current may be written as 

= 5{N^'(x, t)) = J d=^xi J diie"(*i-*) (^7V^(x, t), C(xi, ti)^ . (62) 

Equations (61) and (62) describe the response of the fluid to small perturba- 
tions. 

In the next step we need to match the averages of the conserved quantities 
(61) and (62) to those that appear in the hydrodynamic equations. For 
this purpose, we must first specify the parameters /3 and /i as the inverse 
temperature and the chemical potential in local equilibrium, respectively. In 
order to enforce such matching it is sufficient to require that in the rest frame 
of the fluid [1, 3, 30, 31], 



(f°°(x,t)) 
(7V°(x,t)) 



(T°°(x,t)),, 
(^°(x,i)),. 



(63) 
(64) 
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or in a covariant form u^5{T^^)vy — and u^S{N^) — 0. Secondly, we 
need to identify the parameter with the four-velocity appearing in our 
hydrodynamic equations in Sec. 2. In general, such identification requires a 
specific frame, e.g., the Landau-Lifshitz (LL) frame, where the four- velocity 
is parallel to the energy flow, or the Eckart frame where the four-velocity 
is parallel to the charge flow. From Eqs. (61)- (62) we see that in the frame 
where = 0, neither 6{T^'^) nor 6{N^) necessarily vanishes, i.e., we require a 
transformation to either the LL or Eckart frame. We choose to work in the LL 
frame, in which case the three- velocity should be boosted as = 6{T^^)/[e + 
P), where £ = (f°0),, P = ^i{f'''')i/3, and n = {N^)i are the equilibrium 
values of energy density, pressure, and charge density, respectively. (In the 
Eckart frame the boost is given by = S{N^)/n.) 

Written in the LL frame, the hydrodynamic stress tensor and electric 
current are [3] 

rpi^u ^ {f^''')i + S{f'"') -uVS{fp^)A'"' -u''ufS{fp^)A'"', (65) 
= (7V'^)^ + 5(7VM)-^^M''5(T,,)A'^^ (66) 

where the quantities on the left are those appearing in the hydrodynamic 
discussion of Sec. 2. According to Eqs. (63) and (64) we have u^5(^^y)u^ = 
and u^S{Nfj) — 0, so any linear combination of these expressions can be freely 
added to the right-hand side of Eq. (65). This allows us to rewrite T^^'^ and 
n'^ as 

rpnu = {f^"')i + S{K^'''), (67) 
n^" = {Nt^)i-S{G'') (68) 

after introducing 

+ (e^A^, + $aS^,) , (69) 

G, = j-^UpT'^'^A^.-N'^A,,, (70) 

where 0^, 0^, and $a are defined in Appendix B. 

The advantage of casting T^'^ and in the form given by Eqs. (67) and 
(68) is that the operator C can be decomposed as (see Appendix B) 

C = GpWa + pkp^w'"', (71) 
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where Gp has odd spatial parity and Kp„ has even spatial parity. Later it 
will be seen that such a decomposition permits us to write the transport 
coefficients in a symmetric manner. 

For further convenience, we decompose the energy-momentum tensor T^^ 
and the current TV'' as follows: 

Ni' = fiu^^ji^^ilf^ (72) 
the expansion coefficients being defined by 



e 












Pi. 






P\\ 












Ri" 












ft 












i 







2" " / -'po-' 



(73) 

The operators h^, R^, ir^^'^, and Jf^ are all perpendicular to u^, while Rf^, tt^^'^, 
and J^^ are also perpendicular to b^. The physical meanings of these operators 
are obvious. With these definitions we finally obtain 

= n^'' - Pj_E^'' + P\\b%'' - R^b'' - R'b^', (74) 

= ^hTp^'^ - -^^^ - ^^M' (75) 

with P± = P±- {Q/3 + $/?)£ - (©a + ^a)n, and = P|| - e^e - e„n. 

The dissipative terms can now be written as correlation functions con- 
structed from Eqs. (74) and (75), 



J J —00 
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X (K^,{yi,t),KP^{Ki,ti)^Wp^{xi,h), (76) 
/(x,t) = 5(iV^(x,t)) 

^ -J d^xi J diie^(*i-*) (G"^(x,i),G'^(xi,ii)) V,a(xi,ii), 

(77) 

where we have invoked Curie's principle that the correlators between opera- 
tors with different spatial parity must vanish. 

Now suppose that the changes of the thermodynamical forces are suffi- 
ciently small over the correlation lengths of these expressions, such that they 
can be factored out from the integral. Then we obtain the desired linear 
relations between the thermodynamical forces and the irreversible currents. 
The coefficients in these linear relations are easily identified as the viscosities 
and thermal conductivities: 

r)"""'' = pjd'^xij dhe'^''-'^ (^K^''{x,t),kf"'{xi,h)y (78) 
k"'' = -pjd^^ij dhe'^*'-^^ (d^'{x,t),G''{xi,ti)y (79) 

For practical purposes, it is more convenient to express these transport coef- 
ficients in terms of retarded Green's functions, which can be computed using 
the methods of equilibrium Green's functions. Straightforward manipula- 
tions lead us (see Appendix Appendix C for details) to 

^H^P<r ^ ^^G'^-Pa^o,^^) , (80) 



du! 



a)->0 

(81) 



where G^^'^" {}<i,uj) and G^'^(k,a') are the Fourier transforms of the retarded 
Green's functions 

Gr(x,t) = -i^(t)[G'^(x,t), 6^(0,0)]. (82) 

(Note that we are using the same notation for Green's functions in real space- 
time and for their Fourier transforms in momentum-energy space.) Next, we 
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match Eq. (80) with Eq. (37), which leads us to 

A^,V'""'%A = -3C||, 

A^.T^'^'^^'^A,^ = 10,^0 + + 47^2 + 3C± + 3C||, 

bMabpv'^'"'^ = ^?7o + ?7i + 3C||, 

bf^ab.b^r'"''' = 27^4. (83) 

These equations uniquely determine the viscosity coefficients; explicitly, they 
are given by 

a = ^A^.jy'^-^S^^, 
Cii = -iA^,<-^6,6^, 

7)0 = ^(-A^,A„^-2A^,M/3 + 2A^„A,;3+4A^aM/3 + MA&/3X""'^, 
o 

771 = - (A^,A„^ + SA^Mfs - 2A^^A,0 - AA^Khp + bb^b^bM r]^'"'^ , 
m ^ \ (A^,A«^ + 2A^,6„6^ - 2A^«A,^ - SA^J),hp - hh^h,hj)p) 7]^^""^, 

V4 = lb,ab,b0v'"'"f'. (84) 

We obtain the thermal conductivities in similar fashion by matching Eq. (81) 
with Eq. (38). Their exphcit form is given by 



2 



«x = -i^b^aK^r (85) 
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Finally, we wish to express the transport coefficients in terms of the Green's 
functions. Substituting the results for the viscosity and thermal conductivity 
tensors, Eqs. (80) and (81), into Eqs. (84) and (85), we arrive at 



a 

Cii 

^0 

Vi 

V3 



2Gl B(0,a;) + G'|B,(0, 



ld_ r 

3 duj 

1 d r 

3 dio 
1 d 

--— ImGf ^^.(0,cu) 



0,a;) + G| . (0,c^) 



8 du '^'^^ ' ' 

1 d 



0,a;) 



(86) 



where the retarded Green's function is defined as 



G\^{^,t) = -iQ{t) i(x,i), 5(0,0) 



In deriving Eq. (86), we have used the properties — ImG'^^(w, = ImG^ 
and G^^{u,b'^) = €a€bG^^{u, —6°") of the retarded Green's function, which 

are valid for any Hermitian bosonic operators A and B. (Here (e^) is 
the parity of the operator A (B) under time- reversal.) These relations are 
the consequence of the Onsager's reciprocal principle [1]. Similarly, for the 
thermal conductivities we have 



d 

1 d 

2 OLO "/"'t^'^.t' uj^o 



0)^0 



(87) 



For practical purposes it is sufficient to compute the transport coefficients in 
the rest frame of the fluid. In this case the correlation functions for the 77 's 
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and k's simplify, and we find 




(88) 



Equations (84)- (88) are our main results. They express the transport coeffi- 
cients of strongly magnetized fluids in terms of equilibrium (retarded) corre- 
lation functions, which can be computed by methods of equilibrium statistical 
mechanics. They take into account the fact that the magnetic field breaks 
the translational symmetry of the problem, thereby increasing the number of 
transport coefficients compared to the nonmagnetic case. These expressions 
are valid at large field strengths (but below values at which vacuum quan- 
tum fiuctuations become important) and can be used to compute transport 
coefficients in quantizing fields. A specific example was given in our earlier 
work [22]. 

It is instructive to compare our results to those for an isotropic, un- 
magnetized {B = 0) fluid. In the isotropic case, the most general tensor 
decompositions for the viscosity tensor and thermal conductivity tensor are, 
respectively. 




(89) 



The corresponding Kubo formulas are 



c 



d_ 
du! 



lmGfp{0,uj) 
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^ = l^^^G^..G..M (90) 



3 dco G"'^" 
or when written in the rest frame, 



C = -— ImG|^(0,a;) 



d_ 

K = -AiniGg^^^^(o,a;) , (91) 



?7 = -^ImG'|i2^i2(0,a;) 



a)->o' 



where 



T^"" = ( A A - A I fP'^ 



\de J ^ \dn J 
P = (92) 



These expressions are obtained directly from Eq. (86) and Eq. (88) by taking 
the hmit 5 = 0. However, it is to be noted that even in the isotropic case 
one has {P±,P\\) ^ {P,P) for pressure-pressure correlation functions. 

Thus, we have succeeded in expressing the transport coefficients of a 
strongly magnetized fluid in terms of certain correlation functions of quan- 
tities appearing in the energy-momentum tensor. For practical evaluation 
of these correlation functions, one needs to specify a suitable Lagrangian 
describing the system and a many-body scheme which allows for resum- 
mation of polarization-type diagrams in a controlled manner. Typically, 
non-perturbative resummation schemes are needed to obtain correct infrared 
behavior of the transport coefficients. Concrete applications have been pre- 
sented, for example, in Refs. [32, 33]. 

4. Summary and Discussion 

In this work we have extended the development of covariant magneto- 
hydrodynamics of relativistic fluids initiated in Ref. [22]. A hydrodynamical 
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description of physical processes in highly magnetized compact objects can be 
achieved starting from conservation laws for conserved charges and a suitable 
expansion scheme that systematically resolves the large-length-scale and low- 
frequency limits. In the most general case, the dissipative function contains 
five shear viscosities, two bulk viscosities, and three thermal conductivity 
coefficients, which encode the transport processes at the microscopic level. 
By utilizing Zubarev's method of non- equilibrium statistical operators, we 
have related these transport coefficients to correlation functions of quantities 
entering the energy-momentum tensor in equilibrium. This is done to linear 
order in the expansion of non-equilibrium statistical operators with respect to 
the gradients of relevant statistical parameters (temperature, chemical poten- 
tials, and momentum.) The linear perturbation theory recovers the ordinary 
dissipative Navier-Stokes theory. In general, in order to obtain stable and 
causal hydrodynamics one needs to expand to second order with respect to 
the gradients, in which case additional transport coefficients arise. It would 
be an interesting task to compute these coefficients in the magnetohydrody- 
namic theory in the limit of large magnetic fields. 

We have shown how to relate the correlation function, constructed from 
ingredients of the energy-momentum tensor, to the retarded Green's func- 
tions of equilibrium theory at vanishing momentum. The static coefficients 
are obtained, as usual, by taking the limit of vanishing frequency. In order to 
obtain an explicit form for the transport coefficients, the energy-momentum 
tensor (or, equivalently, the Lagrangian) describing the system needs to be 
specified. The two-particle response functions computed at any finite order in 
the perturbation theory do not give the correct dependence of the transport 
coefficients on the couphng(s) of the underlying theory. Therefore, a specific 
resummation scheme is needed in each case. The Kubo formulas obtained 
in this work provide the necessary starting point for implementing such a 
program for various systems affected by strong magnetic fields. 
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Appendix A. Dependence of 6'^"6^^ + 6'^'^6^" on combinations in 
Eq. (32) 

In a previous paper [22], we have treated the combination h^'^lfl^ ^ift^-l^jji^a 
as an independent tensor structure, since it has the symmetries of the vis- 
cous tensor rj^""'^. Here we show that, in fact, it can be written as a hnear 
combination of the tensors appearing in Eq. (32). By direct calculation we 
find 



9iJ-v 9nP 9tip 9a<j 

gau gap gap g au 

g\u gxp gxp gxa 

g^v g(i3 g^p g^a 

= 9ij.u9a0 - gtiisgua + bi3 {b^g^iu - bf,gi,a) - Up {Uagp.u - u^gua) 

-bu {bagp./3 - bp.gal3) + {Uagp.13 " M^S'a/?) 
- {bij.Ua - baU^) {b^Ufs - bfjU^) 
— ^nu^aP — ^ua^nl3 + babp^nu + b^b^Aa^ — bijbpi!i>.i,a — bi^ba^p.^- 

(A.1) 



Hence we have 
bp^bv^ + bi^abnp 



= 2A^,A„^ - A,«A^^ - A^„A,^ + 26„6^A^, + 26^6,A„^ 

—bpb/sA^a — biyba^np — bvbpAfjoi — bnbaA„p 
= 2(i)-(ii) + 2(iii)-(v). (A.2) 



Appendix B. Derivation of Eq. (71) 

The purpose of this Appendix is to give the details of the derivation of 
Eq. (71). Wc first cstabUsh the following lemmas. 

Lemma 1: The derivatives of /3 = 1/T and a = fijj, in thermal equilibrium 
can be written as 



Da = -P{ej + ^a(l>): 



(B.l) 
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where 



,dP\ ^ ^fdM 



n,B 



Proof: First, using Pd/j, — da — jid^, dP — —P'^dT, and e + P — Ts + fin we 
have 

dP = sdT + ndn + MdB 

= -T{e + P)dl3 + Tnda + MdB, (B.3) 

or, equivalently, 

fl^l ._^(, + P), fff) fl^l ^M. (B.4) 



Secondly, using ds — Pde — adn + ^MdB, we obtain 



dl3\ _ _ /da 



\ fm ^ [ d{PM) \ 



Based on the above thermodynamic relations, we may now write 

- = (a,/- (n).- 



e.n 



d(3dP d[5dP\ 



de d/3 dn da J 



B 
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_ rdP\ „ ^ fdM 



From these equations, we immediately recover Eqs. (B.l). 
Lemma 2: If the electric field is neglected, we have 

Wd^a = i±L{TKDu^ -}fd,T), (B.8) 
ni ^ 

^'""d^a = i±^{T^P-Du,-^'"'d,T), (B.9) 



or equivalently. 



= i±Z (TD«^ _ v^T) . (B.IO) 



Proof: When the electric field is neglected, — 0, we have <9^T^m = 
-n^F^"' = nE'' = 0. Substituting T^"^ = ^^^(m^m'^ - S^^ - h^'h'') for the 
energy-momentum tensor of the electromagnetic field, a direct calculation 
shows that 

S^^ (a^ In S - Du"" + IfdpU") = 0. (B.ll) 
From the first Maxwell equation, dy [B^^u^ — B^u^) = 0, we have 

u^dy {Bi'u'' - B^u^) = = u^'n'^d^B^ - d^B^ , (B.12) 
which implies 

dyU" = -U'dy In S - hi'Du^. (B.13) 
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Contracting the second Maxwell equation, d^H'^'^ — n^, with b^", we obtain 

= b'"'buud^H + HbP''df'b, 



= E^'^d^H + HE^^ib'^d^b^-Du^), (B.14) 
which implies 

E,p}fd^¥ = E^pDu^ - E^^d'' In H. (B. 15) 
Substituting Eq. (B.ll) into Eq. (B.15), we obtain 

E^pb^d^bP = E,p {DuP - dP In M) . (B.16) 
The equation of motion of an ideal fluid is 

^ — 'JixJ- — ^M-' FO 

= u^De + (£ + P^)u''e + (£ + P^)Du'' 

-E'^^'d^P^ + MBib^d^b^ + + b^^d^P^. 



(B.17) 



Contracting this equation with b^ and Spj, yields 



{e + Pi_%Du'' - MBdpb'' - b^^d^P^ = 0, (B.18) 
{e + P±)Ep,Du'' - Ep^d^'P^ + MBEp.b^d^b'' = 0. (B.19) 

Substitution of Eq. (B.13) into Eq. (B.18) followed by some manipulation 
gives 

{e + PA.)b^Dvy + MBilfd^ In B + b"" Du^) - b^'i.sdpT + n9^/x + Md^B) = 0. 
We then have 

nb'^d^^i = (/xn + TsyDu^ - sb'^d^T, (B.20) 

or equivalently, 

b-'d^a = {Tb^Du" - b-'d^T) , (B.21) 
27 



thereupon proving Eq. (B.8). Substituting Eq. (B.16) into Eq. (B.19), we 
obtain 

or equivalently, 



which proves Eq. (B.9). 

Making use of Lemmas 1 and 2, we may now write 

= /3u^A^pDuP + 13V ^u^ + u^u^DI3 + u^A^pd^ji 
+ P {A^pKa + Qpu„u^Ap^ + ^pUp^u^Ep^) d^u", 

-I3up {Q^Ap, + d'u". 
Substituting these relations into 

C = f^,{^,t)d^/3''{^,t)-N^{^,t)d^a{^,t), 
we immediately obtain 

C ^ (u^f'^'A.p - ii^iV^A,^^ {TDuP - VT) 
+ P T^'^iA^pA^^ + e^Ui^u^Ap^ + ^pu^u^Ep^) 



(B.22) 
(B.23) 



+N''Up {Qa^pa + ^^a^pa) 

n 



e + P 



'■vp 



+ /3 



f^" {AppAya + QcjUpU^Apa + ^j3UpUyEp„) 



+N''up {e^Ap^ + $«S 



per J 



pa 



with wf" = (V'u'" + VuP)/2. 



(B.24) 

(B.25) 
(B.26) 



(B.27) 
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Appendix C. Integration of Kubo Correlator and Reteirded Green's 
Function 

In this Appendix, we derive the relation between the space-time integral 

of the Kubo correlator and the retarded Green's function. This relation was 
used to derive Eqs. (80)- (81). Let us denote the integral of a general Kubo 
correlator by 



1 = Jd^^'J* diV(*'-*)(x(x,t),F(x',t'))- 



(C.l) 



First, we rewrite the Kubo correlator in the following form, assuming that 
limt'^-oo{X{x,t)Y{x',t'))i = {X{x,t))i{Y{x',t'))i so that the correlation 
vanishes at — )■ — oo: 



X(x,t),F(x',t') 



rfr(X(x, t) e-^"F(x', t')e''^ - (F(x', t')) 



dr 



ds 



1 
1 

^(X(x, t) [y (x', s + ir)- (y(x', s + ir))^ )i 



ds{ X(x,i),f(x',s) 



(C.2) 



In order to obtain the last equality we have used the Kubo-Martin-Schwinger 
relation {X{t)Y{t' + = {Y{t')X{t))i. 
Then we find that 



I = 



/ dt'e'^''-'^^ / ds{\x{x,t),Y{^',s) 



- I r dt'e'^''-'^^ f dsGR{-K-^,t- s) 

J J —oo P J —oo 

r dt'e'''\ f ds f — e*"MimG'i^(k,cu) 

-oo P J-oo J 27r k->0 



duj 



1 



27r iuj(iuj + e) k->o 



lim GR{k,uj) 
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T- 7 7 1™ lim 

p J 2tt iuj[iuj + e) w->o k-*-o 



- ^ lim lim —^^(k, a;), (C.3) 

the contour integration being performed in the upper-half plane where the 
retarded Green's function is analytic. \{ A = the real (imaginary) part of 
the retarded Green's function is an even (odd) function of cj, 

7 = -4 lim lim ^^^^(k, uj). (C.4) 

This result is used in Sec. 3 to obtain Eqs. (80)- (81). 
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